16 



Table 1. Samples from the IOC source catalogue, which is available in full online at www.mrao.cam.ac.uk/surveys/10C. The extraction procedures for the various parameters are described fully 
in Paper I. The columns are the source name; the group name, which indicates overlapping sources (the group is named after the source with the highest peak flux density in the group and the 
number of objects in the group is given in brackets); the right ascension, Opk, and declination, <5p]j, (J2000) of the source peak; the peak flux density, Spk, and its uncertainty, SSp]^; the right 
ascension, aju, and declination, Sj^, of the source centroid (J2000; relevant for extended sources); the integrated flux density. Sin, a-nd its uncertainty, SSi^; the critical value, e^riti of the major 
axis, above which the source is defined as extended, having deconvolved the telescope's point-source response from the fitted 2-D Gaussian; the lengths, emaj a-nd e^i^, of the major and minor axes 
and position angle, eg, of the source after the deconvolution - the emin and eg columns are left blank for sources classified as point- like; the source type, t - point-like (P) or extended (E); a flag, 
'*', indicating that the error from the approximation of the point-source response is significant and, consequently, the results from the Gaussian fitting ought to be treated with caution; and the 
region of the survey in which the source falls - deep (D), shallow (S) or neither (N). Note that the parameters from the Gaussian fitting are omitted if the fitting did not converge. 
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ABSTRACT 

In a previous paper, the observational, mapping and source-extraction techniques used for 
the Tenth Cambridge (IOC) Survey of Radio Sources were described. Here, the first resuhs 
from the survey, carried out using the Arcminute Microkelvin Imager Large Array (LA) at an 
observing frequency of 15.7 GHz, are presented. The survey fields cover an area of » 27 deg^ 
to a flux-density completeness of 1 mJy. Results for some deeper areas, covering » 12 deg^, 
wholly contained within the total areas and complete to 0.5 mJy, are also presented. The 
completeness for both areas is estimated to be at least 93 per cent. The IOC survey is the 
deepest radio survey of any significant extent (> 0.2 deg^) above 1.4 GHz. 

The IOC source catalogue contains 1897 entries and is available online from the survey 
website (www.mrao.cam.ac.uk/surveys/10C). The source catalogue has been combined with 
that of the Ninth Cambridge Survey to calculate the 15.7-GHz source counts. A broken power 
law is found to provide a good parameterisation of the differential count between 0.5 mJy and 
1 Jy. The measured source count has been compared to that predicted bv lde Zotti et alj ( l2005h 
- the model is found to display good agreement with the data at the highest flux densities. 
However, over the entire flux-density range of the measured count (0.5 mJy to 1 Jy), the 
model is found to under-predict the integrated count by s; 30 per cent. 

Entries from the source catalogue have been matched to those contained in the catalogues 
of the NRAO-VLA Sky Survey and the Faint Images of the Sky at Twenty Centimetres survey 
(both of which have observing frequencies of 1 .4 GHz). This matching provides evidence for a 
shift in the typical 1.4-to-15.7-GHz spectral index of the 15.7-GHz-selected source population 
with decreasing flux density towards sub-mJy levels - the spectra tend to become less steep. 

Automated methods for detecting extended sources, developed in the earlier IOC paper, 
have been applied to the data; x 5 per cent of the sources are found to be extended relative to 
the LA synthesised beam of ^ 30 arcsec. Investigations using higher-resolution data showed 
that most of the genuinely extended sources at 15.7 GHz are classical doubles, although some 
nearby galaxies and twin-jet sources were also identified. 



* We request that any reference to this paper cites 'AMI Consortium: 

Davies et al. 2010' 

t Email: m.davies@mrao.cam.ac.uk 

t Email: t.franzen@mrao.cam.ac.uk 



1 INTRODUCTION 
I.I Background 

The Ninth Cambridge ( 9C) Survey of Radio Sources 
dWaldram et al] l2003l l201(]h . carried out using the Ryle Tele- 
scope (RT) at an observing frequency of 15.2 GHz, was a 
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milestone in the exploration of tlie higii-radio-frequency sky, as 
the first survey of significant extent and depth at such a high radio 
frequency. Since the publication of the first 9C paper, extensive 
survey work has been carried out using the Australia Telescope 
Com pact Array at 20 GHz (ATCA; Ricci et al. 2004; Sadler et al.' 
12006; Massardi et al. 2008, 2010; Murphy et al. 2010). The two 
surveys are complementary, with 9C probing deeper flux-density 
levels (down to 5.5 mjy) and the ATCA surveys covering shallower 
and wider areas (most recently, the whole southern sky). 

It is well known that high-frequency radio surveys are highly 
time-consuming. The scaling of interferometer primary-beam areas 
with frequency (oc v"^), and the typical synchrotron spectra of radio 
sources (oc v""-^) conspire so that the time required to carry out a 
survey of equivalent depth and sky-coverage, using a telescope of 
fixed aperture diameter, scales as v^ *. Things are somewhat better 
if it is assumed that the available bandwidth scales linearly with fre- 
quency. However, the fact that the noise temperatures of the avail- 
able front-end, low-noise amplifiers used in interferometers tends 
to increase with frequency must also be taken into account. 

For these reasons, relatively little survey work has been un- 
dertaken at high radio frequencies and the knowledge of the source 
population remains poor. Nevertheless, familiarity with the prop- 
erties of this population is important for the interpretation of the 
results from observations of the Cosmic Microwave Backgr ound 
(CMB), such as those made by Planck fauber et a"l]|2010l) . At 
mm wavelengths, foreground radio sources are the do minant source 
of co n tamination of small-scale C MB anisotropies l lde Zotti et al.l 
ll999l) . IWaldrametalJj2003l l20ld) have demonstrated that extrap- 
olation of the flux densities of sources at low frequencies cannot 
be relied upon to predict their high-frequency properties, which 
emphasises the value of survey work at the frequencies of interest 
(> 10 GHz) for CMB work. 

Samples of bright sources selected at high radio frequen- 
cies have sig nificant proportions with flat or rising spectra (see, 
for e xample, iTavlor et alj I2OOII ; I AMI Consortium: Davies et al.l 
I2OO9I ). In the main, these sources are believed to be blazars 
with synchrotron self-absorbed spectra; the self-absorbed compo- 
nents of such sources are often highly variable (see, for example, 
IaMI Consortium: Franzen et al.ll2009() . 

High-frequency-selected samples also include appreciable 
numbers of sources with convex spectra, peaking at GHz fre- 
quencies (see, for example, Bolton et al. 2004). Some of these 
GHz peaked spectrum (GPS) sources (see lO'Deall 19981 for a re- 
view) are believed to be associated with young objects, which 
later expand into powerful radio sources, though many are domi- 
nated by^emission from a strongly-beamed self-absorbed compo- 
nent IBolton et al. 2006). Surveys such as the 9C provide flux- 
density-limited samples, which are useful for gaining further un- 
derstanding of the evolution of such objects. 



1.2 This work 

Since the 9C survey was carried out, the RT has been trans- 
formed, by the installation of new front-end receivers and 
back-end electronics (including a new correlator), into the Ar- 
cminute Microkelvin Imager (AM I) Large Array (LA) (see 
I AMI Consortium: Zwart et alj|2008l . for a detailed description of 
the telescope). The LA is a radio synthesis-telescope, located a; 
19 m above sea level near Cambridge. It is used to observe at 
a centre frequency of 15.7 GHz and has a usable bandwidth of 
4.5 GHz. At this frequency, the telescope has a full-width-at-half- 



maximum (FWHM) primary beam of a; 5.5 arcmin and a resolution 
of a 30 arcsec. 

The LA has been used to carry out the Tenth Cambridge (IOC) 
Survey of Radio Sources. As part of this survey, the improved flux 
sensitivity of the LA, compared with the RT, has been used to ex- 
plore the 15-GHz-band sky to sub-mjy levels. In a previous paper 
( IAMI Consortium: Franzen et alJl2010L hereafter Paper I) detailed 
technical information regarding the survey strategy, mapping and 
source-extraction techniques for the IOC survey was provided. In 
this paper, the first results from 10 fields, including the 15.7-GHz 
source count, are presented. Throughout this paper any equatorial 
coordinates use equinox J2000 and spectral indices are defined us- 
ing the convention that S oc v^" . 



2 THE IOC SOURCE CATALOGUE 

The techniques used for observing, mapping and source extraction 
are described fully in Paper I. The fields, the positions of which 
are given in Section [3] were surveyed using a 'rastering' tech- 
nique, with observations being carried out between 2008 August 
and 2010 June. Each field was observed using a set of telescope 
pointing directions spaced at 4.0 arcmin intervals and lying on a 
2-D hexagonally-gridded lattice, projected on to the plane of the 
sky. A raster map that combines the individual CLEANed maps be- 
longing to each of the pointing directions was produced for each 
field; the raster map for one of the survey fields is shown in Fig.[T] 
In addition, a noise map that shows how the noise varies across the 
raster map was created for each field; these noise maps are used in 
identifying sources, as described in Paper I. 

Information about the sources was extracted from the raster 
maps, using a combination of in-house software and tasks belong- 
ing to the Aip£], for inclusion in the IOC source catalogue. Source- 
finding was carried out using a flux-density threshold of 4.62cr; the 
reason for this slightly unusual choice is explained in Sections [2. 11 
andlOl 

In Section [3] areas within each of the survey fields, bounded 
by lines of constant right ascension and declination, complete to 
flux densities of I.O and 0.5 mJy, are defined. A short section of 
the catalogue is shown in Table[T] The methods used to extract the 
various parameters are described fully in Paper 1 ; the final column 
indicates which of the areas (see Section[3]( each of the sources lies 
within. The complete source list, which contains 1897 entries, is 
available online at www.mrao.cam.ac.uk/surveys/10C. It is eventu- 
ally planned to also make the fits maps for each of the survey fields 
publicly available on this website. 

Individual positional error estimates have not been assigned to 
each source. The positional errors in both RA and Dec. for a source 
detected at the 5-cr level, in any of the survey fields, are estimated 
to be 3-4 arcsec; this range reflects the fact that the synthesised 
beam is slightly elliptical and has dimensions which vary with field 
declination. The assessment of the source positional accuracy was 
made using simulations in which point sources were inserted into 
the map of one of the survey fields; the extracted source positions 
were compared with the nominal values. 

The results of the simulation were found to agree well with 
the positional errors that would be expected from theory, taking ac- 
count of Gaussian thermal noise. Higher-resolution follow-up ob- 
servations are required to assess the positional accuracy for the 

' ASTRONOMICAL IMAGE PROCESSING SYSTEM - www.aips.nrao.edu/ 
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Figure 1. The raster map of one of the IOC survey fields (J0024+3152). The map was produced by combining approximately 1000 individual constituent 
maps, weighted according to their RMS noises; the individual maps were CLEANed before being combined. The area of lower noise in the centre of the raster 
map is clearly visible. 



brighter sources. Only the highest-flux-density source in the IOC 
survey catalogue has a counterp art in the catalogue of the Jodrell- 
VLA Astrometic Survey (JV AS ; IPatnaiket al]| 19921 : iBrowne et al] 
ll998l : IWilkin"son et aljl99^ . The measured positions of this source 
agree within 3.5 arc sec. 

There were some problems with the analysis and mapping of 
the IOC survey data that required some adaptations to the source- 
finding procedure described in Paper I. These adaptations are de- 
scribed in the remainder of this section. 



2.1 Checking the flux-density scale 

As a check of the raster maps' flux-density scale, a; 50 of the bright- 
est sources detected in the maps were selected for pointed follow- 
up observations, carried out during 2010 June and August using the 
LA. To avoid complicating the analysis, only sources that showed 
no evidence of extension in the raster maps were selected. The 
data from these observations were mapped using the same CLEANing 
scheme used for the raster maps but phase self-calibration was also 
applied. 

Fig. [2] shows the peak flux densities for each of the sources 
measured using the raster and (self-calibrated) pointed maps. The 
figure indicates that the flux densities of the sources measured 



from the raster maps are systematically low compared with those 
from the self-calibrated pointed observations. However, because 
the pointed follow-up observations were carried out about 2 yr af- 
ter the commencement of the raster observations, it is important to 
consider the efl'ect of source variability on this result. 

Fig. [2] indicates that the pointed flux density of one source is 
almost twice its raster flux density; this diff'erence is almost cer- 
tainly attributable to genuine flux-density variability. Having said 
this, the number of genuinely variable sources within the sample 
is likely to be small. The sources selected for pointed follow up 
have flux den sities ranging between approximately 10 and 40 mjy. 
Results from lWaldram et al.1 (|2010') indicate that 15-GHz-selected 
samples containing sources with flux densities in this range are 
likely to be dominated by steep-spectrum sources, which do not 
ordinarily display significant variability. 

Nevertheless, the median percentage difference has been used 
to quantify the discrepancy between the pointed and raster flux den- 
sities because it is less sensitive to genuine source variability than 
the mean, the value of which could be strongly affected by a small 
number of highly-variable sources within the sample. The median 
percentage difl'erence between the pointed and raster flux densities 
was calculated to be 8.2 per cent with an uncertainty in this value 
of » 2 per cent. 
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on timescales much longer than 10 min. If the measured phase for 
a calibrator source does change by more than 30° between succes- 
sive visits to the source, the affected data are automatically omitted 
during the data-reduction stage; however, such a large phase dif- 
ference between successive calibrator visits is observed only very 
rarely. Similarly, any data for which the estimated error in the mea- 
sured phase of the calibrator is > 15° are also omitted automati- 
cally. 

The discrepancy between the flux densities measured using 
the phase-self-calibrated pointed observations and the raster obser- 
vations is instead attributed to phase errors in AMI data resulting 
from the u neven spacings of t he time delays in the telescope's lag 
correlator. iHolIer et alj | |2007|) explain this problem in detail and 
propose a solution. In practice, however, further work is required to 
analyse the data from the correlator correctly. 



^ ' 1 1 1 1 ' ' 1 

5 10 15 20 25 30 35 40 

Peak raster flux density/mJy 

Figure 2. The peak flux densities of » 50 sources as measured from the 
raster maps and from maps created using self-calibrated data from pointed 
follow-up observations. Points lying on the line represent sources having 
identical flux densities when measured using the raster and pointed maps. 



Since a number of the IOC survey fields overlap with areas 
mapped as part of the 9C survey, as an additional check the flux 
densities measured from the pointed observations were compared 
with the values measured from pointed 9C observations. Having 
corrected the phase errors in the IOC pointed observations using 
self-calibration, there was found to be good agreement between the 
9C and IOC pointed values with a median percentage difference of 
< 1 per cent - the median is again used for the same reason as 
given above. It is noted that, owing to the small difference in the 
observing frequencies and the typical spectra of radio sources, the 
IOC flux densities might be expected to be slightly lower than the 
9C values - this was, in fact, the case. 

Data from a large number of observations, carried out during 
a range of weather conditions, were combined to create the raster 
maps; whereas the data used to produce the pointed maps were 
collected during single, short observations during relatively good 
(dry) observing conditions. As a result, the weather conditions dur- 
ing the pointed observations may not be typical of those during the 
raster observations. However, in practice data collected during pe- 
riods of poor weather are either omitted entirely or significantly 
downweighted with respect to data collected during good weather 
conditions. 

A modulated noise signal, injected at the front end of each an- 
tenna, was monitored throughout the raster observations. The data 
were then weighted, on a sample-by-sample basis, according to the 
value of the 'rain gauge' - that is the ratio of the power of the mod- 
ulated noise signal to the total power input to the correlator (which 
is kept constant) to that obtained in cool, dry, clear weather con- 
ditions. In addition, data for which the rain gauge was less than 
50 per cent (this is rather a conservative criterion) were omitted. 
This is because the amplitude correction applied to the data be- 
comes unreliable during heavy rain, as no attempt to account for 
atmospheric absorption is made in applying this correction. 

Atmospheric-related phase effects are not thought to make 
any significant contribution to the phase errors present in the data. 
During LA observations a phase-calibrator source is observed for 
I min at a time at 10 min intervals. Even using the telescope's 
longest (IlO-m) baseline, the measured phase varies only slowly 



2.2 Correcting the sources' flux densities for phase errors 

Unfortunately, the great majority of the sources in the IOC raster 
maps are detected with insufficient signal-to-noise to allow self- 
calibration to be successfully applied to the data. Therefore, a cor- 
rection to the source flux densities based on the difference be- 
tween the flux densities measured from the raster maps and the 
self-calibrated pointed maps is applied; as a final step in the source- 
finding procedure, all source flux densities are multiplied by 1.082 
before inclusion in the IOC catalogue. 

The uncertainties in the flux densities are increased to take ac- 
count of this scaling and the uncertainty in the correction factor 
(a; 2 per cent), which is in any case small compared with the esti- 
mated LA calibration uncertainty of 5 per cent. As in Paper I, the 
uncertainty in an uncorrected peak flux density, 5, is taken to be 
^jcr^ + (0.055)-, where cr„ is the thermal noise at the source posi- 
tion, estimated from the noise map. The uncertainty in the corrected 
value is, therefore, 

/0-2 + (0.055)2 " 
1.0825 y ° ^+0.022. (I) 

Initially, it was intended to carry out source finding at 5cr; 
since, assuming Gaussian statistics for the map noise, such a 
scheme would result in a highly reliable catalogue with a; 0. 1 false 
detections. Further, the completeness of such a catalogue would 
have been very high (> 94 per cent) at 0.5 mjy and 1 mjy in the 
deep and shallow regions respectively. However, owing to the effect 
of phase errors, a source of 5 = I mJy falling within the shallow 
area of the survey, which ought to be detected at > 5cr, will only be 
detected at > 5cr/ 1.082 = 4.62cr. 

Therefore, in order to achieve the desired high levels of com- 
pleteness at 0.5 mJy and I mJy, it was decided to carry out the 
source finding at A.&lcr. The catalogue completeness is discussed 
fully in Section|4] The slight relaxation in the source-finding crite- 
rion is likely to have only a small adverse effect on the reliability of 
the catalogue. Assuming Gaussian statistics for the map noise, and 
given the number of synthesised beam areas in the survey maps, it 
is estimated that about one source in the final catalogue will be a 
false positive. 

2.3 Excluding areas around bright sources 

The IOC raster maps often display an increased level of noise 
around bright (> 15 mJy) sources. This is attributable to amplitude, 
phase and deconvolution errors in the data. The elevated noise level 
close to bright sources is generally not fully reflected by the noise 
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maps, because the noise is higlily non-Gaussian in these regions. 
Therefore, such detections close to sources of 5 > 15 mjy were au- 
tomatically excluded from the final source catalogue. In addition, 
detections close to a number of fainter sources (the faintest being 
x 9 mJy) were also excluded manually from the final catalogue 
after inspecting the maps. 

Empirically, the area of elevated noise around each bright 
source (the 'exclusion zone'), from which the source detections 
were rejected, was found to be well-represented by a circle, cen- 
tred on the source position, of radius 

r. = 12 arcmin, (2) 

where Sp^ is the peak flux density of the source. Table|2]shows the 
centre coordinates and radii of the exclusion zones applied to the 
survey data. 

The dynamic range of the LA - that is the ratio of a source's 
peak flux density to the flux density of the brightest artefact in a 
map close to the source - is « 50:1. The brightest source detected 
in the IOC survey is a 270 mJy. A conservative approach, which 
also serves to simplify the completeness areas for construction of 
the source count, has been taken by assuming that sources with 
S > 9 mJy can be detected anywhere in the total areas. In contrast, 
any putative source with a peak flux density less than this value that 
falls within the exclusion zone of a bright source is not included in 
the final IOC source catalogue. Thus, the area used for calculating 
the source counts for sources with 5 < 9 mJy does not include 
the exclusion zones around bright sources. The total area excluded 
from around bright sources is 0.6 deg^. 



3 SURVEY FIELDS 

The IOC survey fields are centred at J0024-1-3I52, J0259+2610, 
J0734-I-5432, J0824-H6931, J0939H-3I15, J1046H-5904, 
J1052-I-5730, J1524+432I, J1543+4420 and J1733H-4I48. 
Fig. [3] shows the positions of these fields, which were chosen to 
be widely spread in RA and away from the Galactic plane. For 
each field, areas complete to I.O and 0.5 mJy (apart from the 
exclusion zones) have been defined by selecting those areas in 
which the noise, (r„, estimated using the relevant noise map is 
O.I < CTn < 0.2 mJy and cr„ < 0. 1 mJy respectively. Over a large 
portion of the areas, the estimated noise is significantly lower than 
the upper bounds and, therefore, the completeness is close to 100 
per cent for both regions at their nominal completeness levels. 

The areas complete to 0.5 mJy are entirely contained within 
the areas complete to 1.0 mJy and are referred to as the 'deep' re- 
gions. Areas of higher noise, complete to 1.0 mJy, but excluding 
the deep regions, are referred to as the 'shallow' areas. Fig.|4]shows 
the deep and shallow regions for one of the fields. The survey cat- 
alogue includes a flag for each source indicating whether it falls in 
the deep (D) or shallow (S) areas. Sources that fall outside these 
regions altogether (i.e. in areas of higher noise towards the edges 
of the raster maps) are indicated by 'N'. For a source with evidence 
of extension the flag is based on the source's centroid position, oth- 
erwise the peak position is used. The 'total' areas are those areas 
complete to 1 mJy but not excluding the deep regions - in other 
words, the combined deep and shallow areas. The lines of right as- 
cension and declination bounding the total and deep areas, for each 
of the fields, are given in Tables[3]and|4]respectively. 



Table 2. The centre positions of the exclusion zones around bright sources 
and their radii. 



RA 


Dec. 


Radius (arcmin) 


00:20:50.4 


+31:52:29 


3.39 


00:21:29.8 


+32:26:60 


3.63 


00:23:09.9 


+31:14:01 


4.26 


00:26:06.2 


+32:08:33 


2.83 


00:28:10.7 


+31:03:46 


3.85 


00:29:20.4 


+32:16:55 


3.17 


00:29:33.1 


+32:44:58 


4.45 


02:59:55.1 


+26:27:26 


2.31 


03:01:05.5 


+25:47:16 


2.92 


03:01:37.3 


+25:41:54 


3.04 


07:31:17.4 


+53:38:58 


4.25 


07:36:52.9 


+54:29:17 


2.70 


08:18:16.1 


+69:16:53 


4.06 


08:23:02.5 


+69:14:20 


2.94 


09:35:59.5 


+31:27:27 


3.17 


09:36:36.9 


+32:03:35 


3.31 


09:37:06.2 


+32:06:58 


5.41 


09:41:03.2 


+31:26:14 


3.30 


09:41:46.2 


+31:55:03 


3.05 


09:42:08.8 


+32:06:42 


3.48 


10:47:19.3 


+58:21:14 


4.93 


10:49:40.0 


+58:35:31 


3.35 


10:50:07.1 


+56:53:37 


3.38 


10:50:54.0 


+58:32:33 


3.50 


10:51:41.4 


+59:13:08 


3.80 


10:52:25.4 


+57:55:08 


3.46 


10:52:54.5 


+59:22:18 


3.56 


10:54:26.9 


+57:36:48 


3.69 


15:20:41.6 


+44:13:18 


3.43 


15:21:49.4 


+43:36:37 


12.51 


15:27:51.8 


+43:52:05 


2.85 


15:28:19.8 


+42:33:35 


4.01 


15:40:33.5 


+44:34:01 


3.32 


15:41:10.0 


+44:56:34 


4.58 


15:42:23.1 


+43:59:15 


3.81 


15:46:04.5 


+44:49:14 


3.02 


17:25:34.5 


+41:53:03 


4.33 


17:27:49.3 


+42:21:40 


4.45 


17:29:01.9 


+41:40:04 


2.68 


17:30:41.6 


+41:02:58 


6.14 


17:31:23.7 


+41:01:38 


3.01 


17:37:59.6 


+41:54:51 


5.00 


17:38:35.4 


+42:21:43 


3.08 


17:40:08.9 


+41:36:09 


6.31 


17:40:17.2 


+42:14:30 


2.93 


17:40:52.1 


+42:34:47 


5.79 



4 COMPLETENESS 
4.1 Simulations 

Simulations were carried out to investigate the completeness of 
the survey. A number of realisations were used (twelve for the 
deep areas and thirteen for the shallow) in which 250 equal-flux- 
density, simulated point sources were inserted into the raster map 
of J0024+3152. The flux density of the simulated sources was dif- 
ferent for each of the realisations. The positions of the simulated 
sources were chosen randomly but were not altered between re- 
alisations; to avoid the simulated sources affecting each other, it 
was insisted that no simulated source could lie within 2 arcmin of 
any other. Sources were not inserted in the exclusion zones around 
bright sources, since, as explained above, the completeness limit is 
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Table 3. The areas complete to 1 mJy. 



Field 


RA 


range 


Dec 


range 


Area (deg^) 


J0024+3152 


00:19:54.2 


to 


UU.Zy . JO. J 


+31:00:04 




+32"43'05 


3.56 


J0259+2610 


02:56:43.6 


to 




+25:19:17 


to 


+27"02" 17 


2.24 


J0734+5432 


07:30:02.8 


to 


U/.jo.jU. / 


+53:41:42 


to 


+JJ.ZJ.U / 


2.16 


TnS94-(-fiQ3 1 


0R-17-0R 6 


to 


08:31:23.8 


+68-41 -04 


to 


+70:22:41 


2 11 


J0939+3115 


09:36:22.1 


to 


09:42:25.9 


+30:24:22 


to 


+32:06:18 


2.20 


J1046+5904 


10:39:25.5 


to 


10:52:31.2 


+58:13:12 


to 


+59:55:07 


2.86 


J1052+5730 


10:47:36.8 


to 


10:57:21.3 


+56:48:56 


to 


+58:10:22 


1.78 


J 1524+4321 


15:20:10.7 


to 


15:29:24.6 


+42:30:20 


to 


+44:11:21 


2.83 


J1543+4420 


15:38:30.5 


to 


15:47:46.6 


+43:28:43 


to 


+45:11:13 


2.83 


J1733+4148 


17:25:37.0 


to 


17:41:01.3 


+40:57:34 


to 


+42:40:06 


4.90 



Table 4. The areas complete to 0.5 mJy. Note that there are three deep regions associated 
with J0259+2610. 



Field 


RA 


range 


Dec 


range 


Area (deg'^) 


J0024+3152 


00:20:59.1 


to 


00:28:27.1 


+31:19:38 


to 


+32:23:29 


1.69 


J0259+2610 


02:57:50.5 


to 


03:01:24.1 


+25:38:31 


to 


+26:43:51 


0.87 




02:56:50.2 


to 


02:57:50.5 


+25:53:05 


to 


+26:28:22 


0.13 




03:01:24.1 


to 


03:02:31.4 


+25:53:05 


to 


+26:28:22 


0.15 


J0734+5432 


07:31:47.0 


to 


07:37:06.6 


+54:00:59 


to 


+55:04:59 


0.82 


J0824+6931 


08:19:21.5 


to 


08:29:26.1 


+68:57:47 


to 


+70:07:19 


1.02 


J0939+3115 


09:37:32.8 


to 


09:41:15.9 


+30:42:52 


to 


+31:47:13 


0.85 


J1046+5904 


10:41:30.0 


to 


10:50:24.5 


+58:32:47 


to 


+59:35:08 


1.19 


J1052+5730 


10:49:57.7 


to 


10:55:32.7 


+57:08:03 


to 


+57:51:28 


0.54 


J 1524+4321 


15:21:32.3 


to 


15:27:57.3 


+42:50:56 


to 


+43:52:21 


1.19 


J1543+4420 


15:39:51.6 


to 


15:46:21.3 


+43:49:49 


to 


+44:50:48 


1.18 


J1733+4148 


17:26:44.0 


to 


17:39:49.8 


+41:17:01 


to 


+42:19:40 


2.54 




32°30'- 



Figure 3. The positions of the IOC survey fields, shown using an equatorial- 
plane projection with the North Pole at the centre. The declination circles 
are at intervals of 10° and the Galactic plane is indicated by the thick black 
line. 



expected to be much higher for these areas compared with the re- 
mainder of the maps. In order to investigate the completeness as a 
function of flux density, the flux density of the sources was changed 
between each realisation. 

For each realisation, the ordinary source-finding procedures 
were applied to the map and the proportions of simulated sources 
that were recovered was calculated for the shallow and deep ar- 
eas. Three of the simulated sources were found to lie too close to 
real sources to be detected separately. In these cases, the simulated 




31°00 



0''30™ 28 



Figure 4. The shallow (lighter shading) and deep areas belonging to 
J0024+5432. 



source was considered to be detected if the recovered source posi- 
tion was closer to the simulated rather than the real source position. 

Fig. |5] shows the proportion of sources detected for each of 
the realisations compared with the expected detection rate, which 
is calculated for a particular region of the map as follows. For a 
specific source flux density, the probability of detection can be cal- 
culated at the position of each pixel by taking account of the value 
of the noise map at that position and by using Gaussian statistics. 
The probability of detection over the whole region is straightfor- 
wardly calculated by averaging over pixels. 

The plots show that, in general, there is good agreement be- 
tween the predicted completeness curves and the results of the 
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Figure 5. Results of a simulation to investigate the completeness of the 
survey. Filled circles show the proportion of the simulated sources recov- 
ered as a function of flux density within the deep (cTn ^0.1 mJy) area 
of J0024+3152. The solid line shows the completeness predicted based on 
the noise-map pixel values assuming Gaussian statistics for the noise. The 
open squares and dashed line show the results for the shallow (0. 1 < (Tn < 
0.2 inJy) area. 



simulation. However, for the fainter sources, the detection rate 
is slightly higher than predicted, whilst the converse is true for 
the brighter sources. Source confusion is likely to be responsible 
for these results. A number of the highest flux-density simulated 
sources, that would otherwise have been detected, were not recov- 
ered owing to their proximity to real, bright sources. This effect 
prevented the completeness from reaching 100 per cent as quickly 
as predicted. At fainter flux densities, however, the opposite effect, 
whereby a source is boosted in flux density owing to its proximity 
to a real, fainter source, so that it is unexpectedly detected, becomes 
important. 

The peak flux densities of sources detected with low signal-to- 
noise ratios (SNRs) are typically biased slightly high, because the 
peak position tends to be coincident with a positive noise fluctua- 
tion. This is an additional factor serving to boost the detection rate 
at low SNRs. The presence of this effect was confirmed by extract- 
ing the flux densities of the sources at the precise positions with 
which they were simulated; the values were found to be systemati- 
cally low compared with the extracted peak flux densities and were 
found to reflect better the simulated flux densities. 



4.2 Using the noise maps to estimate the completeness 

Having established by simulation that, assuming Gaussian statistics 
for the noise, the noise maps can be used to provide reasonable 
estimates of the survey completeness, the noise maps from all the 
fields were used to estimate the completeness of the IOC survey. 
The probability of a source of true flux density 5 being included in 
the survey catalogue when located on a pixel with a corresponding 
noise-map value of (T„ was taken, according to Gaussian statistics, 



0.2 0.4 0.6 0.8 1 1.2 1.4 

Source flux density/mjy 

Figure 6. The estimated probabihty of detection for the shallow (dashed 
line) and deep (solid line) areas for all survey fields. The dot-dashed hori- 
zontal line indicates a probability of one. The dot-dashed vertical lines indi- 
cate the 0.5-mJy and 1.0-mJy nominal completeness limits for the deep and 
shallow areas respectively. 



to be 



P{S > 4.62o-„) : 



J4.62o-„ 



(x- 5 /1. 082) 
exp — dx. 



(3) 



In carrying out this calculation the fact that sources are detected 
with flux densities lower than their true values has been taken into 
account by including the factor of 1/1.082. 

Knowing the actual distribution of noise-map pixel values. 
Equation [3] can be used to estimate the completeness, as a func- 
tion of flux density, for the shallow and deep regions of the survey; 
Fig.|6]shows these estimates for both areas. For the shallow regions, 
the catalogue is estimated to be » 94-per-cent complete by 1 mJy 
and 99-per-cent complete by a 1.16 mJy. For the deep areas, the 
catalogue is estimated to be » 93-per-cent complete by 0.5 mJy 
and 99-per-cent complete by « 0.61 mJy. 



5 SOURCE COUNTS 

The 15.7-GHz differential source counts have been calculated by 
binning the sources from the final catalogue according to their peak 
flux densities, except for those that display evidence of being ex- 
tended relative to the telescope synthesised beam (this is the case 
for 5.5 per cent of sources), for which integrated flux densities were 
used. 

The binned differential-source-count data are given in Table[5] 
For the highest flux-density bin, data from the entirety of the to- 
tal areas have been used. The bins for sources with flux densi- 
ties between 1 and 9 mJy include sources from the total regions 
but exclude the areas around bright sources given in Table |2] The 
bins for sources with flux densities between 0.5 and 1 mJy include 
sources from the deep regions (again excluding areas around bright 
sources). The source count is not calculated for 5 > 25 mJy, since 
it is biased low in this flux-density range; using 9C data, several of 
the fields were selected to contain as few sources with 5 > 25 mJy 
as possible. 
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Table 5. Data for the IOC source count. 



Bin start 


Bin end 


Number of 


Area 


5 (mJy) 


S (mJy) 


sources 


(deg2) 


9.000 


25.000 


46 


27.46 


5.500 


9.000 


51 


26.86 


2.900 


5.500 


142 


26.86 


2.050 


2.900 


135 


26.86 


1.500 


2.050 


148 


26.86 


1.250 


1.500 


113 


26.86 


1.000 


1.250 


160 


26.86 


0.900 


1.000 


36 


11.96 


0.775 


0.900 


56 


11.96 


0.680 


0.775 


51 


11.96 


0.600 


0.680 


64 


11.96 


0.540 


0.600 


61 


11.96 


0.500 


0.540 


46 


11.96 



At 0.5 mJy, the completeness limit of the deep areas, the sur- 
vey is limited by thermal rather than confusion noise. Above this 
flux-density level, it is estimated that there are typically 170 LA 
synthesised beam areas per source. 

The effect of the calibration errors of a 5 per cent is negligible 
serving to boost the number of sources in each bin by << 1 per cent 
(owing to the sign of the slope of the source count). However, bias 
teddingtorj 1 1 9 1 3h due to the thermal noise will play a more im- 
portant role for the bins at the faintest flux-density levels. Given 
the slope of the counts and the noise properties of the deep region, 
the number of sources in the faintest flux-density bin is expected to 
be boosted by a 7 per cent. However, this effect is almost exactly 
balanced by the small degree of incompleteness that affects these 
faintest flux-density bins. Consequently, no corrections for the ef- 
fect of Eddington bias or incompleteness have been applied. 

The IOC differential source count is shown in Fig. [7] As in 
all subsequent plots showing source counts, Poisson errors, on 
the number of sources in each bin, are indicated for each of the 
points and the bars parallel to the flux-density axis represent the 
bin widths, not error bars. 

An attempt was made to fit a single power law to the data. 
However, such a model did not appear to fit the data well. Conse- 
quently, a broken power law has been fitted to the data. A method 
was used whereby the sum of the squared differences between 
the measured and predicted areas under the curve, over all bins, 
was minimised. In carrying out the minimisation, the points were 
weighted according to their respective Poisson errors. The positions 
of the points in the S direction within the bins have been plotted, in 
Fig.|7] on the basis of the fitted exponents to reflect the 'centre of 
gravity' of each bin. The fitted differential source count is 



n(S) = 



dN 
dS 




Jy-' sr-' for 2.8 < 5 < 25 mJy 
Jy-' sr-' for 0.5 < 5 < 2.8 mJy. 



Tests were carried out to check that the broken power-law 
model does indeed provide an improved fit to the data, compared 
with the single power-law model. An F-test indicated that the null 
hypothesis - that the data are more likely to have been drawn from 
the simpler model - can be rejected at > 99.9 per cent confidence. A 
comparison of the models using Akaike's information criterion was 
similarly emphatic, indicating that the broken power-law model is 
~ 1000 times more likely to be correct than the simpler model. 

In order to check the self-consistency of the source count, 
the data were used to construct separate counts for the shallow 
and deep regions. Fig. [8] shows the counts from the two regions 



1e+08 



100000 



27.46 deg; ^ 
26.86 deg, " 
11.96 deg^ - 



0.001 



0.01 



S/Jy 



Figure 7. The IOC differential source count. Poisson eiTors, on the number 
of sources in each bin, are indicated for each of the points. The bars parallel 
to the flux-density axis represent the bin widths, not error bars. The fitted 
broken-power-law count is indicated by the solid line. 



1e+09 



1e+07 



1e+06 



Deep " -<5 - 
Shallow ^ -■- - 




S/Jy 

Figure 8. Differential source counts from the shallow and deep regions of 
the IOC survey. The fitted IOC source count is indicated by the solid line. 



overlaid. The deep counts are derived from an area amounting to 
1 1.96 deg- (except for the highest flux-density bin, which contains 
data from 12.19 deg^). The shallow counts use data from 14.90 deg^ 
(15.27 deg^ for the highest flux-density point). The plot shows good 
agreement, within the uncertainties, between the counts derived 
from the two regions over the common flux-density range. 

5.1 Adding in data from the 9C survey 

It is possible to extend the IOC source count to higher flux densi- 
ties by the inclusion of data from the 9C survey. It is also possible 
to improve the source count statistics between 5.5 and 25 mJy by 
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Table 6. Data for the combined 9C and IOC source counts. 



0.8 



Bin start 
S (mjy) 


Bin end 
S (mJy) 


Number of 
sources 


Area 
(deg2) 


500.000 


1000.000 


8 


520.00 


200.000 


500.000 


27 


520.00 


100.000 


200.000 


47 


520.00 


60.000 


100.000 


92 


520.00 


40.000 


60.000 


97 


520.00 


30.000 


40.000 


99 


520.00 


25.000 


30.000 


79 


520.00 


16.000 


25.000 


62 


124.60 


12.000 


16.000 


64 


124.60 


10.000 


12.000 


48 


124.60 


9.000 


10.000 


15 


47.83 


6.400 


9.000 


48 


47.23 


5.500 


6.400 


44 


47.23 


2.900 


5.500 


142 


26.86 


2.050 


2.900 


135 


26.86 


1.500 


2.050 


148 


26.86 


1.200 


1.500 


140 


26.86 


1.000 


1.200 


133 


26.86 


0.900 


1.000 


36 


11.96 


0.775 


0.900 


56 


11.96 


0.680 


0.775 


51 


11.96 


0.600 


0.680 


64 


11.96 


0.540 


0.600 


61 


11.96 


0.500 


0.540 


46 


11.96 



including 9C data, since the 9C survey contains regions that are 
complete over this flux-density range. There is a small difference 
in the observing frequencies of the 9C (15.2 GHz) and IOC surveys 
(15.7 GHz). Therefore, in combining the data sets, the source flux 
densities from the 9C survey catalogue have been corrected to take 
into account this difference. The correction was made by assuming 
a typical spectral index between 15.2 and 15.7 GHz that varies as a 
function of source flux density. 

The assumed flux-density-dependent spectral index is indi- 
cated in Fig. [9] The correction was calculated by fitting a logarith- 
mic function to the points in the plot. For each of the three lowest 
flux-density points (indicated by filled squares), the median flux 
density versus the median value of the 1.4-to-15.2-GHz spectral 
index, □'[4", of sources, with flux densities in the relevant ranges, 
detected as part of the 9C survey (see Table 9 of Waldram et a l. 
I2OI0I) has been plotted. For the highest flux-density point (indi- 
cated by the filled triangle) , the median flux density versus median 
value of Q-|g, measured by IAMI Consortium: Davies et al.l ( |2009|) . 
for sources belonging to a flux-density-limited source sample 
jLopez-Carriego et al.ll2007D from the 3-yr Wilkinson Microwave 
Anisotropy Probe {WMAP) data has been plotted. The number of 
sources belonging to the samples represented by these points ranges 
between 84 and 381. 

For both the 9C and WMAP samples it has simply been as- 
sumed that the spectral indices can be extended over the small ad- 
ditional frequency ranges to 15.7 GHz (from 15.2 GHz in the case 
of the 9C samples and from 16 GHz for the WMAP sample). The 
approach followed is admittedly not perfect, particularly since the 
sample for the highest flux-density point was selected at 33 GHz. 
Nevertheless, because the difference in observing frequencies be- 
tween the surveys is very small, the corrections are similarly small 
(at most a few percent) and the method is considered acceptable. 

Fig. 1101 shows the combined 9C and IOC di fferential source 
count. Data from areas of the 9C s urvey presented in lWaldram et al.l 
1 I2OO3I) and lWaldram et al.l ( I2OI0I) have been used in constructing 
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Figure 9. The filled squares represent the median 1.4-to-15.2-GHz spectral 
indices for complete 15.2-GHz-selected source samples from the 9C survey. 
The filled triangle represents the median 16-to-33-GHz spectral index for a 
complete 33-GHz-selected sample from the WMAP 3-yr data. The dashed 
line indicates the assumed typical spectral index, as a function of flux den- 
sity, used to make coirections to the flux densities of individual 9C sources, 
to account for the small diflerence in the observing frequencies of the 9C 
and IOC surveys. This allowed the IOC source count to be extended to flux 
densities > 25 mJy by the inclusion of 9C data. The typical spectral index 
was calculated by fitting a logarithmic function to the data points. 



the count. The data used for each bin are shown in Table [6] For 
the intermediate flux-density ranges, for which there are data from 
both surveys, some of the 9C survey data were excluded. This was 
to avoid double counting the areas that were surveyed as part of 
both the 9C and IOC surveys. 

Using the same method as described above, a broken power 
law was fitted to the binned differential count. As previously, a bro- 
ken power law was found to be a significantly better fit to the data 
than a single power law. Again, an F-test indicated that the null hy- 
pothesis - that the data are more likely to have been drawn from the 
simpler model - can be rejected at > 99.9 per cent confidence. The 
best-fit broken-power-law parameterisation of the source count is 

_ dA' (48 (f )^ " Jy"' sr-' for 2.2 mJy < 5 < 1 Jy 
d5 I 340 (^) Jy-' sr-' for 0.5 < 5 < 2.2 mJy. 

This fitted count is found to give a good fit to the data, with a re- 
duced chi-squared value of 0.75. The probability of obtaining a 
reduced chi-squared value greater than 0.75 by chance, given the 
number (20) of degrees of freedom, is 78 per cent. The fit is indi- 
cated in Fig. 1101 

5.2 Comparison with the de Zotti model 

In Fig. [TT] the combined 9C and IOC source count is compared 
with the latest ver sion of the 15-GHz source-count model by 
Ide Zotti et al.l ( l2005h . extracted from their websit^l on 201 1 March 
01; for completeness, the model counts, over the relevant flux- 
density range, are provided in Appendix |A] No attempt has been 

^ http://web.oapd.inaf.it/rstools/srccnt/srccnt_tables 
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Figure 10. The combined 9C and IOC 15.7-GHz differential source count. Dif ferent symbols are use d to distinguish between the areas t hat were used to derive 
the count for the various flux-density bins; '9Cr indicates areas presented in lWaldram et alj i2003h . '9C2' refers to areas presented in lWaldram et alj j201(ll) 
and 'IOC is used to designate areas presented in this paper The fitted count is indicated by the solid line. 



made to correct for the small frequency difference between the 
measured and model source counts but this is likely to make lit- 
tle difference to the overall conclusions. 

The model count is in good agreement with the measured 
count at the high-flux-density end. However, the shape of the plot- 
ted model count is somewhat different from that of the measured 
count; so that, with decreasing flux density, the model first over- 
predicts and then, below a; 5 mjy, under-predicts the measured 
count. 

The total number of sources per steradian with flux densities 
between 0.5 mJy and 1 Jy, predicted by the model, was calculated 
by integrating the model difi'erential source count between these 
limiting flux densities. Since the predicted counts are given at a 
number of discrete flux densities (see Table IaT] ). the integration 
was carried out piecewise by approximating the model count as a 
power law between successive pairs of points. 

The number of sources per unit area predicted by the model 



was found to be only 70 per cent of the measured value. Because the 
differential count is largest at the low flux-density-end, the under- 
prediction of the count at the lowest flux densities dominates over 
the over-prediction at slightly higher flux densities, explaining the 
30-per-cent deficit over the entire range. 



6 MATCHING WITH 1.4-GHz SURVEYS 

The complete and unbiased sample of sources from the deep areas 
of the survey (i.e. those sources with flux densities between 0.5 and 
25 mJy in these areas) was matched with entries from the 1.4-GHz 
catalogues of the Faint Ima ges of the Radio Sky at Twenty Cen- 
timetres (FIRST; iBecker. Wh ite, . ^Hgjf and 1995) su rvey and the 
NRAO-VLA Sky Survey (NVSS; ICondon et al.l 199^ . The sample 
comprises a; 650 sources in total. There are NVSS data for all the 
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Figure 11. The normalised (S^ **^) combined 9C and IOC differential 
source count. Symbols indicating the areas from which the counts were 
derived are identical to those in Fig. |10| The fitted broken-power-law pa- 
rameterisation is indicated by the so lid line. The dashed line indicates the 
prediction of the latest version of the lde Zotti etsd\ i2005l) model. 



IOC survey fields; however, some of the IOC survey fields are not 
covered or are only partially covered by FIRST. 

The NVSS, which has a resolution of 45 arcsec (slightly lower 
than that of the IOC survey) has a completeness of 50 per cent at 
~ 2.5 mjy and 99 per cent at x 3.4 mjy. The FIRST survey has 
higher resolution of 5 arcsec and is complete to 1 mJy. Owing to 
the significant difference in resolution between the IOC and FIRST 
surveys, the flux densities of any FIRST sources lying within one 
LA synthesised beam of a IOC source were summed for compari- 
son with the IOC flux density. 

Some results from matching with the NVSS and FIRST cata- 
logues are shown in Tables|7]and[8]respectively. In both cases, the 
sample has been divided into various 15.7-GHz flux-density bins 
and the percentage of sources with flux densities falling in each 
bin without a corresponding match in the relevant low-frequency 
catalogue is indicated. A limiting 1.4-to-I5.7-GHz spectral index, 
ttiim, based on the completeness limit of the relevant low-frequency 
catalogues (3.4 mJy for NVSS and 1.0 mJy for FIRST), is also in- 
dicated for each bin. Any source with a spectral index greater than 
this limiting value ought to be detected in the low-frequency cata- 
logue. 

In addition, the percentage of sources with a > 0.5 and a > 
0.81 is given for each bin. For the matching to the NVSS catalogue, 
for all but the three highest flux-density bins, the percentages of 
sources with a > 0.5 ought to be regarded as lower limits. This is 
because the values of q-u,^, for the lower flux-density bins, are > 0.5. 
The choice of the second value of 0.81 avoids this problem, since 
for all flux-density bins and for matching to both low-frequency 
catalogues, all sources with O' > 0.81 should be matched. 

It is noted that the percentages of sources with o- > 0.5 and 
a > 0.81 from matching to the FIRST survey are systematically 
lower than those from matching to the NVSS. This is most likely 
due to the fact that extended emission detected as part of the NVSS 
is resolved out as part of the FIRST survey. 



An interesting efi^ect, whereby the source population shifts to- 
wards a flatter-spectrum population with decreasing flux density, is 
observed in the data. This is most clearly seen in the matching to 
the NVSS from the change in the fraction of sources with a > 0.81. 
As mentioned previously, the fraction of sources with a > 0.5 rep- 
resents a lower limit for a number of bins and so the evolution in 
this fraction cannot readily be interpreted as evidence for a change 
in the population. In the matching to the FIRST survey, the effect 
is best demonstrated by the change in the fraction of sources with 
a > 0.5. Owing to the fact that FIRST coverage is not available 
for all the IOC fields, the statistics for the fractions of sources with 
a > 0.81 are poorer compared with those for the matching with the 
NVSS. Also, the interpretation of the results is complicated by the 
mismatch in resolution noted above. 

By matching so urces detected as part of the 9C survey to 
the NVSS catalogue, IWaldram et al.l ( I201Q) observed a shift in 
the I5-GHz-band source population from being dominated by flat- 
spectrum to being dominated by steep-spectrum sources with de- 
creasing flux density between » 100 mJy and « 10 mJy - this 
efl^ ect is illustrate d by Fig. [9l A similar effect has been observed 
by iMassardi et aU t201Q) by matching entries from their 20-GHz 
source catalogue to low-frequency catalogues. Their large sample 
size allows them to trace the spectral evolution of the source pop- 
ulation in detail between 40 mJy and 1 Jy. The median value of 
af^ is found to become rapidly larger with decreasing flux density 
below K 80 mJy (E. Sadler; private communication). 

The IOC data suggest that, at fainter flux densities, this trend 
is reversed with a move back towards a flatter-spectrum popula- 
tion. A shift towards a flatter-spectrum population has been ob- 
served in_sub;mJ^_data_^om low-frequency surveys (see, for ex- 
ample, |Prandonrernj200^ and references therein). 



7 EXTENDED SOURCES 

Using the methods described in Paper I, 4.8 per cent of the IOC 
sources are classified as extended and 5.9 per cent as overlapping 
(the categories are mutually exclusive). Some of the 'extended' 
sources might be classified as such owing to source confusion or ef- 
fects arising from noise. Similarly, chance juxtapositions of sources 
are likely to result in sources with no genuine association being 
classified as overlapping. Consequently, simulations in which sim- 
ulated point sources were inserted, at randomly-selected positions, 
into the IOC survey fields were carried out in an attempt to pro- 
vide an estimate of the fraction of IOC sources that are genuinely 
extended relative to the LA synthesised beam. 

In total « 1600 simulated point sources were added into the 
IOC data in the uv plane. The sources were simulated with flux 
densities drawn randomly from the measured source count. As for 
the simulation to assess the survey completeness, no simulated 
source was allowed to lie within 2 arcmin of any other. The stan- 
dard mapping and source extraction procedures, described in Pa- 
per I, were applied to the data. Of the simulated point sources that 
were detected, 1.4 ± 0.5 per cent were classified as extended and 
2.2 + 0.4 per cent were categorised as overlapping, where the uncer- 
tainties have been estimated from the scatter in the results between 
the different survey fields. 

The results from the simulation suggest that x 100 x (4.8 - 
I.4)/4.8 « 71 per cent of the sources classified as such are gen- 
uinely extended and that a 100 x (5.9 - 2.2)/5.9 ~ 63 per cent 
of overlapping sources are actually a component of a larger ob- 
ject. In the great majority of cases, overlapping sources are found 
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Table 7. Some statistics from matching witfi the NVSS source catalogue. Numbers and percentages of sources with a > 0.5 for the lowest four flux-density 
bins are lower limits. 
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Bin end 
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Unmatched 


a\\m 


Percentage 


a > 0.5 


Percentage 


a > 0.81 


Percentage 


S (mJy) 


S (mJy) 










Unmatched 




a > 0.5 




a > 0.81 


5.0 


25.0 


55 


54 


1 


-0.15 


2 


38 


69 


27 


49 


2.5 


5.0 


65 


62 


3 


0.14 


5 


47 


72 


30 


46 


1.5 


2.5 


99 


87 


12 


0.35 


12 


64 


74 


45 


46 


1.0 


1.5 
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76 


38 


0.52 


33 


>67 


>59 


29 


25 


0.8 


1.0 


79 


50 


29 


0.61 


37 


>42 


> 53 


22 


28 


0.6 


0.8 


125 


60 


65 


0.73 


52 


>59 


> 47 


26 


21 


0.5 


0.6 


107 


43 


64 


0.81 


60 


>43 


> 40 


19 


18 



Table 8. Some statistics from matching with the FIRST source catalogue. 



Bin start Bin end Total Matched Unmatched a\\^^ Percentage a > 0.5 Percentage a > 0.81 Percentage 

S (mJy) S (mJy) Unmatched a > 0.5 a > 0.81 

5l) 251) 36 36 -0.67 17 47 8 22 

2.5 5.0 46 44 2 -0.38 4 23 50 7 15 

1.5 2.5 56 50 6 -0.17 11 30 54 8 14 

1.0 1.5 67 50 17 0.00 25 24 36 10 15 

0.8 1.0 50 41 9 0.09 18 18 36 5 10 

0.6 0.8 88 56 32 0.21 36 35 40 15 17 

0.5 0.6 68 32 36 0.29 53 18 27 7 10 



to have just two components. Therefore, if the two categories are 
taken together - that is, sources consisting of a single compo- 
nent that is extended relative to the telescope synthesised beam 
and extended sources consisting of several overlapping components 
- the fraction of IOC sources that are genuinely extended rela- 
tive to the LA synthesised beam of » 30 arcsec is estimated as 
(5.9 - 2.2)/2 + (4.8 - 1.4) a 5 per cent (the factor of 1/2 re- 
sults from the assumption that the overlapping sources consist of 
two components). This result assumes no clustering of sources on 
small (< 1 arcmin) angular scales. It is noted that a similar frac- 
tio n (x 6 per cent) of ext ended sources was measured from 9C data 
bv lWaldram e7^ ( |2003|) . who carried out the classification on the 
basis of manual inspection of 9C contour maps. 

Contour plots of all the sources identified by the automated 
source-extraction procedures as extended or overlapping were in- 
spected by eye. Sources thought most likely to be genuinely ex- 
tended from their morphologies (for example overlapping sources 
with similar peak flux densities) were identified. IOC contour plots 
of 39 such sources were compared with contour plots from the 
NVSS and, where available, from the FIRST survey. In addition, 
the NASA Extragalactic Databas^ (NED) was searched for objects 
detected within » 2 arcmin of each of the extended IOC sources. 
NED reports entries from catalogues belonging to a number of low- 
frequency radio surveys, in addition to entries from the catalogues 
of optical, infrared and X-ray surveys. 

From these investigations 29 of the 39 sources were found 
to be genuinely extended and two artefacts of source confusion. 
In eight cases it still remained unclear whether the sources were 
genuinely extended. Most of the genuinely extended sources were 
found to be classical doubles but three are nearby galaxies and 
four are twin-jet/tailed sources. One of the twin-jet sources has 
z = 0.267 and another has z = 0.391; the respective centroid po- 
sitions of these sources are 15:41:05.0 -)-43:27:00 and 10:44:56.7 



' http://nedwww.ipac.caltech.edu 
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+59:25:38. Some results for one source of particular interest are 
provided below. 



7.1 8C 0821-1-695 

It was noticed, whilst manually inspecting the contour plots of two 
sources flagged as extended, that these sources were the lobes of 
a 'linear triple', a broad-emission-line radio galaxy or radio-loud 
quasar (see Fig.ll2t. The central component and outer lobes span an 
angular distance of 6.6 arcmin and appear as three separate sources 
in the IOC catalogue, since no linking emission is detected using 
the IOC dat a. The sour ce, 8C 0821-H695, wa s noticed in the 8C 
survey data teeeslll990h bv iLacv et al.l ( Il993h . who measured the 
size of the source to be 1.5 h"' Mpc and to have a redshift of 0.538. 
They estimate the age of the source to be ~ 3x10' yr and argue that 
the source is large because it has expanded rapidly in a low-density 
environment rather than because it is old. Assuming no variability, 
the IOC data impl y spectral steepen ing of the emission from the 
central component; lLacvetalJfl993l) found a spectral index for the 
central component of 0.23 ± 0.04 between 1.4 and 5 GHz. The IOC 
measurements imply a spectral index for the central component of 
0.79 ± 0.06 between 5 and 15.7 GHz. 



8 CONCLUSIONS 

The AMI LA has been used to carry out the IOC survey, the deep- 
est radio source survey of any significant extent (> 0.2 deg^) above 
1.4 GHz. The resulting deep 15.7-GHz source counts are useful for 
the interpretation of CMB data, for which foreground radio sources 
are an important contaminant. The source catalogue also pro- 
vides an invaluable resource for the study of faint high-frequency- 
selected radio sources. The survey covers « 27 deg^ complete to 
1 mJy and a 12 deg- (wholly contained within the larger area) 
complete to 0.5 mJy. The number of sources with 5 > 25 mJy, ap- 
pearing in the survey catalogue, is biased low; several of the survey 
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Figure 12. Contour plots of 8C 0821+695, a giant radio galaxy lying in 
one of the IOC survey fields, from the IOC survey at 15.7 GHz (top) and 
the NVSS at 1.4 GHz. In both cases a factor of V2 separates the contours, 
which start at ±0.1 and ±1.0 mJy beam"' for the IOC and NVSS plots 
respectively. The negative contours are shown using dashed lines. Crosses 
in the IOC plot indicate the position of IOC sources. The cross in the NVSS 
plot gives the position of the central component source at 15.7 GHz. 



fields were chosen to minimise the number of such sources appear- 
ing within the survey areas. In total, 1897 sources appear in the 
4.62-cr catalogue; the faintest being « 100 yuJy. A number of the 
key conclusions from the work are listed below. 

(1) The IOC differential source count was parameterised using a 
broken power law. This was found to provide a significantly better 
fit to the data than a simpler, single-power-law parameterisation. 
The best-fit broken power law is 

AN [ 24 (f)"' ^' Jy-' sr-' for 2.8 < 5 < 25 mJy 
d5 1^ 376 (^) Jy-' sr-' for 0.5 < 5 < 2.8 mJy. 



(2) After having applied corrections to the individual source flux 
densities measured as part of the 9C survey to account for the 
small dilference in frequencies between the 9C and IOC surveys, 
the 9C and IOC data were combined. The addition of the 9C data 
allowed the calculation of the best estimate of the 15.7-GHz dif- 
ferential source count by improving the source-count statistics for 
5.5 < 5 < 25 mJy and by providing data for a complete sample 
of sources with 25 mJy < 5 < 1 Jy. Again, a broken-power-law 
parameterisation was found to ofi'er a significantly improved fit to 
the data compared with that provided by a single power law. The 
best-fit broken power law is 

_dN I 48 (f)^^^ Jy-' sr-' for 2.2 mJy < 5 < 1 Jy 
d5 y 340 (^) Jy-' sr-' for 0.5 < 5 < 2.2 mJy. 

(3) The model counts bv lde Zotti et al. I( l2005h are found to display 
good agreement with the 9C and IOC data at the high flux- 
density-end of the measured count. However, with decreasing flux 
density the model first over-predicts and then, below about 5 mJy, 
under-predicts the measured count. By integrating the model 
differential source count, the model was found to under-predict 
the total number of sources, with flux densities between 0.5 mJy 
and 1 Jy, per unit area by approximately 30 per cent. This deficit, 
over the entire flux-density range, is attributable to the model 
underestimating the count at the lowest flux densities. 

(4) Entries from the IOC source catalogue were matched to those 
contained in the catalogues of the NVSS and FIRST survey (both 
of which have observing frequencies of 1.4 GHz). The matching 
revealed a shift in the typical I.4-to-15.7-GHz spectral index of 
the I5.7-GHz-selected source population with decreasing flux 
density towards sub-mjy levels. When matching to NVSS, 49 per 
cent of sources with 5.0 < Sioc < 25.0 mJy were found to have 
a}" > 0.81. However, for sources with 0.5 < 5ioc < 0.6 mJy 
the corresponding figure is 18 per cent. The observed trend is in 
contrast to that measu red for sources with higher flux densities by 
IWaldram et al] ( |2010|) . They found that the typical spectral index 
became steeper for sources with decreasing flux densities between 
a; 100 mJy and « 10 mJy. A similar effect, to that measured 
using the IOC data, has been observed as part of lower-frequency 
surveys - in lower-frequency-selected samples significant numbers 
of flatter spectrum sources sta rt to enter at sub-mJy levels (see, for 
example, IPrandoni et al.ll20o3 , and references therein). 

(5) Automated techniques for identifying extended sources, de- 
scribed in Paper I, have been applied to the data. The proportion 
of sources that are extended relative to the LA synthesised beam 
of a 30 arcsec is k 5 per cent; this is similar to the proportion 
of » 6 per cent measured by Waldram et al. (2003). A subset of 
39 extended or overlapping sources, thought likely to be genuinely 
extended on the basis of their 15.7-GHz morphologies, were inves- 
tigated further using higher-resolution data. These data confirmed 
that at least 29 of these sources are genuinely extended; most were 
identified as classical doubles but three are nearby galaxies and four 
are twin-jet sources. 
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Table Al. 15-GHz model source counts by Ide Zotti et al] )2005t) , pre- 
sented as logarithmic, normalised differential counts - specifically 
logioCS'/^dW/dS). Predicted counts are given for a number of different 
source populations; the flat-spectrum counts are made up from the sum of 
the flat-spectrum-radio-quasar (FSRQ) and BL-Lacertae (BL Lac) counts, 
and the total counts are found by summing the flat-spectrum and steep- 
spectrum counts. 



logio S 




Source count (Jy'"^ 


sr"') 




(Jy) 


FSRQ 


BL Lac 


Flat 


Steep 


Total 








spectrum 


spectrum 




-4.000 


-1.563 


-1.020 


-0.91 1 


-1.084 


-0.688 


-3.900 


-1.479 


-0.972 


-0.854 


-0.935 


-0.592 


-3.800 


-1.395 


-0.924 


-0.797 


-0.778 


-0.487 


-3.700 


-1.311 


-0.876 


-0.740 


-0.609 


-0.369 


-3.600 


-1.226 


-0.829 


-0.682 


-0.445 


-0.247 


-3.500 


-1.142 


-0.781 


-0.624 


-0.319 


-0.144 


-3.400 


-1.058 


-0.735 


-0.566 


-0.215 


-0.055 


-3.300 


-0.974 


-0.688 


-0.507 


-0.1 19 


0.030 


-3.200 


-0.890 


-0.642 


-0.447 


-0.028 


0.112 


-3.100 


-0.806 


-0.596 


-0.387 


0.059 


0.192 


-3.000 


-0.721 


-0.550 


-0.326 


0.146 


0.272 


-2.900 


-0.637 


-0.505 


-0.265 


0.232 


0.352 


-2.800 


-0.553 


-0.460 


-0.203 


0.316 


0.431 


-2.700 


-0.469 


-0.415 


-0.140 


0.400 


0.510 


-2.600 


-0.385 


-0.371 


-0.077 


0.483 


0.588 


-2.500 


-0.301 


-0.327 


-0.013 


0.564 


0.666 


-2.400 


-0.217 


-0.283 


0.052 


0.643 


0.742 


-2.300 


-0.133 


-0.240 


0.118 


0.721 


0.818 


-2.200 


-0.050 


-0.197 


0.184 


0.795 


0.890 


-2.100 


0.034 


-0.155 


0.251 


0.868 


0.962 


-2.000 


0.118 


-0.113 


0.319 


0.935 


1.029 


-1.900 


0.201 


-0.071 


0.387 


0.998 


1.093 


-1.800 


0.285 


-0.029 


0.456 


1.056 


1.153 


-1.700 


0.368 


0.012 


0.526 


1.103 


1.205 


-1.600 


0.45 1 


0.052 


0.597 


1.140 


1.249 


-1.500 


0.534 


0.092 


0.668 


1.163 


1.284 


-1.400 


0.616 


0.132 


0.739 


1.163 


1.302 


-1.300 


0.698 


0.172 


0.811 


1.137 


1.305 


-1.200 


0.779 


0.211 


0.883 


1.103 


1.308 


-1.100 


0.860 


0.250 


0.955 


1.076 


1.321 


-1.000 


0.939 


0.288 


1.027 


1.052 


1.340 


-0.900 


1.018 


0.326 


1.098 


1.028 


1.366 


-0.800 


1.095 


0.363 


1.169 


1.003 


1.395 


-0.700 


1.170 


0.401 


1.238 


0.977 


1.428 


-0.600 


1.242 


0.438 


1.305 


0.951 


1.464 


-0.500 


1.310 


0.474 


1.369 


0.924 


1.502 


-0.400 


1.373 


0.510 


1.429 


0.898 


1.541 


-0.300 


1.430 


0.546 


1.483 


0.873 


1.579 


-0.200 


1.479 


0.581 


1.531 


0.850 


1.613 


-0.100 


1.517 


0.616 


1.568 


0.827 


1.641 


0.000 


1.543 


0.651 


1.595 


0.805 


1.660 



pared in Section lS^ Here, the model counts are provided only for 
the flux-density range relevant for comparison to 9C and IOC data. 



APPENDIX A: MODEL SOURCE COUNTS 

TablelAll shows the 15-GHz model source counts bv lde Zotti et al.l 
( 1200511 to which the measured 9C and IOC source counts are com- 
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